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We report the observation of quantum dynamical oscillations of ultracold atomic gases in the F 
and D bands of a single-well optical lattice. We are able to control the Bragg reflections at the 
Brillouin zone edge up to the third order. As a result, we can switch the quantum dynamics from 
oscillations across both the F and D bands to oscillations only within the F-band. Our capability 
to observe these remarkable oscillations comes from the innovative non-adiabatic technique which 
allows us to load ultracold atoms efficiently to the G-band of an optical lattice. 
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I. INTRODUCTION 

There has been a lot of effort both experimentally and 
theoretically to study the quantum dynamics of ultracold 
atoms in optical lattices, such as the celebrated Bloch 
oscillations (BOs) [1-6] and the Landau-Zener (LZ) tun¬ 
neling [3, 7-10]. These studies have focused on the lowest 
band as it is hard to load atoms to high bands and then 
control their quantum dynamics experimentally. Re¬ 
cently, people are pushing the boundary and studying 
quantum dynamics involving more than one bands. The 
effort has resulted in the observation of the Bloch-Zener 
oscillations (BZOs) [11-17], where the quantum oscilla¬ 
tions are between two Bloch bands and the crossing be¬ 
tween these two bands is facilitated by the LZ tunneling. 

However, this kind of oscillations within two Bloch 
bands are very difficult to observe in a simple single-well 
optical lattice. In this kind of simple lattices, the band 
gaps are always smaller for higher bands. If atoms can 
tunnel from the S-band to the P-band, they should also 
be able to tunnel from the P-band to the D band. As a 
result, the oscillations would involve uncontrollably many 
bands [18]. To control oscillations within two bands, one 
has to design lattices with more complex constructions, 
which include binary superlattice of optical waveguide 
arrays [14], honeycomb lattice [15], and mini-band struc¬ 
ture [16, 17]. Here in this work we demonstrate experi¬ 
mentally that we can initiate and control quantum oscil¬ 
lations with two high excited bands in a simple single-well 
optical lattice. Our trick is to use the variable external 
force from the harmonic trap instead of the usual linear 
external potential. 

In this work the Bose-Einstein condensate (BEG) is ini¬ 
tially loaded non-adiabatically into the G-band of a one¬ 
dimensional single-well optical lattice. It subsequently 
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tunnels to the E-band and begins oscillations within the 
E and D bands, which are clearly observed in the momen¬ 
tum space. We can control the Bragg reflection between 
momenta Skhp and —ShkL, which is at the EBZ edge 
between the E and D bands, by tuning the optical lat¬ 
tice strength, kp is the wave vector of the laser forming 
the lattice. When the Bragg reflection is weak, quantum 
oscillations crossing the E and D bands are observed. 
When the reflection is strong, quantum oscillations are 
observed in only the E-band. When the Bragg reflection 
is at intermediate strength, we observe the superposition 
of these two types of quantum oscillations. In our re¬ 
sults, the oscillations crossing the two excited bands can 
last up to 58ms, which is much longer than the BZOs ob¬ 
served in Ref. [17]. During the oscillations, the BEG can 
be displaced up to ±100/im (470 lattice sites) in space, 
which is bigger than the spatial displacement observed in 
super-Bloch oscillations [6] . 

II. EXPERIMENT METHOD TO PREPARE 
THE ATOMS IN THE HIGH BANDS 

In our experiment, a pure BEG of about 1.5 x 10^ 
^^Rb atoms is prepared in a hybrid trap which is 
formed by overlapping a single-beam optical dipole trap 
with wave length 1064nm and a quadrapole magnetic 
trap. The resulting harmonic trapping frequencies are 
{uJx,ojy,uJz) = 27rx (28, 55, 65)Hz. After preparing a BEG 
in the harmonic trap, we use the non-adiabatic shortcut 
method [19-22] to load the BEG in a one-dimensional 
optical lattice (along the x direction) into the G-band 
at the quasi-momentum hq = 0. The optical lattice 
Vq cos^ (kpx) is produced by a standing wave created by 
two counter-propagating laser beams with the lattice con¬ 
stant a = A/2 = 426nm and Vo being the lattice depth. 

This coherent loading method [19, 20] includes a se¬ 
ries of designed standing wave pulses shown in Eig. 1(a). 
It allows us to non-adiabatically load a BEG from the 
ground state of the harmonic trap Ifjo) directly into a 
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FIG. 1: (Color online) (a) The pulse sequence for the loading 
to high bands, (b) TOF pictures taken after loading atoms 
into the S and G bands, respectively. 


target state lipa) = Here |n,g) is the eigenstate of 

Bloch bands, n = S,P,D... is the band index. The pulse 
sequence can be optimized so that the final state is nearly 
a Bloch state (the fidelity can be over 98%). This can 
be verified with the proportion of different momentum 
components or the fast oscillations of different momen¬ 
tum states. For the F and G bands, momenta 

have dominate population while -^2hkL and Ofi/cL for the 
S-band, as indicated by the experimental results shown 
in Fig. 1(b). Due to the symmetry of parity, two series 
of pulses Vo cos^ {kLx) and Vb sin^ (kLx) are needed if we 
want to load directly a BEG into the F-band. For conve¬ 
nience, we load atoms into the G-band with Vb cos^ {k^x) 
pulses sequence in the experiment. 

Our loading process is finished in tens of microseconds 
by applying a series of pulsed optical lattices, and this 
rapid generation of the quantum state in higher bands is 
different from other preparation methods [23, 24]. Our 
method allows us to explore rich physics in high Bloch 
bands predicted theoretically [25-31]. 


III. EXPERIMENTAL OBSERVATION 

After the BEG in the G-band is prepared, we hold 
the lattice and harmonic trap for a period of time t, 
and then switch off all potentials to take pictures af¬ 
ter 28ms time of flight (TOE). Three series of experi¬ 
mental absorption images for the optical lattice depth 
Vb = 5E^, 7.5E^, and IbEr are shown, respectively, in 
Eig. 2(a,b,c). Er = li?k\/2m is the recoil energy with 
m being the atomic mass. The time separation between 
neighboring images in the bE^ series is 1ms and 0.5ms 
in the other two series. These series of images demon¬ 
strate clearly three different quantum oscillations; we will 


explain and analyze them later. Eor convenience, an ex¬ 
tended band structure is drawn in Eig. 2(d), where the 
energy gaps between different bands are marked with Ag 
(s= 1,2,3,4,5,6). 

The BEG is initially loaded in the G-band, where 
the atoms mostly populate equally around two momenta 
EAhkL- As atoms with these two momenta are at the cen¬ 
ter of the trapping potential at the beginning, the only 
possible motion for them is to move either to the left or 
the right. Gonsequently, they lose their momenta while 
gaining harmonic potential energy. This corresponds to 
that the BEG makes a quantum tunneling from the G- 
band to the E-band over the tiny band gap at Ai and 
Aq in Eig. 2 (d). It is impossible for the BEG to move 
up along the G-band due to the conservation of energy. 
When the lattice depth is very high so that the energy 
gap at Ai and Aq is large, it is possible for the BEG 
to stay in the G-band for a long time. We have ob¬ 
served atoms maintaining in the G-band for 5ms when 
Vb = 20Er. Eor the lattice strength of interest in our 
experiment, the BEG always tunnels from the G-band to 
the E-band as soon as the initial loading ends. 

Once the BEG is in the E-band, it continues to lose mo¬ 
mentum while gaining harmonic potential energy. This 
corresponds to that the BEG traverses dynamically along 
the E-band from Ai, Aq to A 2 , A 5 in Eig. 2(d). Once ar¬ 
riving at A 2 and A 5 , the atoms face different ensuing 
dynamics depending on the lattice strength. If the lat¬ 
tice strength is small and the Bragg reflection at A 2 and 
A 5 is weak, the BEG will continue its dynamics into the 
D-band by crossing the band gap. After evolving dy¬ 
namically along the entire D-band, the BEG comes to 
the band gap between D and P bands at A 3 and A 4 . 
This band gap is always large for the lattice strength in 
our experiments. As a result, the atoms at A 3 {—2hkL) 
will be Bragg reflected completely to A 4 (2hkL) while the 
atoms at A 4 will be Bragg reflected completely to A 3 . No 
tunneling to P-band occurs. Afterwards the BEG will re¬ 
verse its dynamics by moving up in momentum from A 4 , 
A 3 to A 5 , A 2 . It eventually arrives at Ae, Ai, finishing 
half of an oscillating cycle. These oscillations crossing 
the two Bloch bands (E and D bands) are driven under 
a variable force from the harmonic trap. They are illus¬ 
trated in Eig. 2 (a) for lattice depth Vb = bEr and their 
period is 24ms. 

Note that the BEG moving up along the bands around 
A 4 , A 3 in Eig. 2(d) by gaining momenta is due to the 
fact that most of the atoms are away from the center 
of the trapping potential and feel an accelerating force. 
This is different from the initial stage when the BEG is 
loaded into the G-band, where most of the atoms are at 
the center of the trap and feel very small force for the 
finite size of the BEG. 

When the optical lattice is strong and the gap at A 2 , 
A 5 is large, the Bragg reflection can dominate the dynam¬ 
ics, forbidding the atoms tunnel from the E-band to the 
D-band. Instead, the atoms at A 2 {—ShkL) will trans¬ 
fer completely to A 5 {ShkL) via Bragg reflection while 
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FIG. 2: (Color online) Quantum oscillations of the BEG in high bands of optical lattices. Shown in (a), (b) and (c) are 
experimental results in momentum space with lattice depth Vo = 5Er, Vb = 7.5Er and IbEr, respectively. The oscillations in 
(a) are across both the F and D bands and oscillations in (c) are only within the F-band; the dynamics in (b) are superposition 
of the oscillations in (a) and (c). (d) is the schematic of extended Bloch bands (P, D, F and G) of a one-dimensional optical 
lattice. The green areas (I and III) are for the F-band and the yellow (II) for the D-band. 


the atoms at will also transfer completely to A 2 . In 
this way, the quantum dynamics is confined within the 
F-band. These oscillations only within the F-band are 
observed in our experiment for Vb = and are shown 
Fig. 2(c) with a period of 17ms. 


When the lattice strength is intermediate, the Bragg 
reflection at A 2 , A^ will be partial: one part of the atoms 
will be reflected and undergo oscillations within the F- 
band; the other part of the atoms will tunnel to the D- 
band and oscillate across both the F and D bands. As a 
result, we should be able to observe a superposition of the 
two kinds of oscillations: across both the F and D bands 
and only within the F-band, when the lattice is at an 
intermediate strength. This is indeed what we observed 
in experiment as shown in Fig. 2(b) and simulated in 
theory as shown in Fig. 3(c) and (d) for Vb = 7.bEr, 
where the two kinds of oscillations are clearly seen, and 
the ratio of them can be tuned by the lattice strength. 


IV. THE THEORETICAL EXPLANATION AND 
SIMULATION 

We have simulated the experiment with the one¬ 
dimensional Gross-Pitaevskii equation (GPE). In 
the simulation, the initial state is 7 /; (x, t = 0) = 

2 pg {x) (j)q=o{x)^ where i^g{x) is the ground state of the 
BEG in a harmonic trap and (j)q=o{x) is the Bloch wave 
function in the G-band at g = 0. We have plotted 
our numerical results in both the real space and the 
momentum space in Fig. 3 for Vq = and Vb = 7.bEr. 
The results in the momentum space agree well with 
the experimental results in Fig. 2. We can infer from 
Fig. 3(b) that the BEG can move away from the center 
of the trap by up to ±100/im for the oscillations in both 
the F and D bands. For oscillations in F-band, our 
numerical results indicate that this displacement can 
be up to ±75/im, which is much larger than what was 
observed in super-Bloch oscillations in Ref. [6]. There 
is a small discrepancy between the theoretical results 
and the experimental results. For example, there is a 
slight off-set in the experimental oscillations around 
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FIG. 3: (Color online) The theoretical results of quantum 
dynamical oscillations for the BEG in high bands in both real 
space and momentum space, respectively. Shown in (a) and 
(b) are results for Vb = 5Er while (c) and (d) are Vb = 7.5Er. 
The white circles in (a) are the result of the semi-classical 
model. 


6-7ms for Vb = bEr] this slight off-set is not seen in the 
corresponding theoretical results. The explanation for 
this small discrepancy becomes clear when we discuss 
Brag reflection in the next section. 

It is difficult to observe oscillations in the real space 
in experiment. To do that, we need to keep the TOF 
very short. Within such a short TOF, the atomic cloud 
is still very dense and the TOF image is not proportional 
to the cloud density. At the same time, the TOF image 
is obscured by the thermal cloud that has no time to 
disperse. 

As the BEG density profile varies smoothly over hun¬ 
dreds of lattice sites, its dynamics can be well described 
by the semi-classical dynamics of Bloch particles [32] , 


, dr 

’'a 

= VqEn{q) 

(1) 

nii 

dt 

= fir) 

(2) 


where En{q) is the nth energy band and f{r) is the 
force acting on the Bloch particle. For Bloch electrons 
in traditional condensed matter physics, we usually have 
f{r) = —eE — x B. In our case, f{r) = —muj‘^xex 
with the harmonic trap frequency 

For simplicity, we describe the energy bands with a 
cosine function as En{q) = ^ cos (qn/kL), with 

n = F, D representing the F and D bands, respectively. 


and \Bn\ the width of the energy bands. For this sim¬ 
plified case, we can solve Eqs.(l,2) analytically and find 
that the oscillation periods are: 

Tp = -^ (3) 

y/mUJx7T 

and 


^ 4hkL f Cl ^ C2 \ 
ymWajTT V-BbJ ’ 


where Ci = f ^ — 7 = 


d(gn) 


= and C 2 = uK{u). 

/cos{qj^a)—cos{qa) 

qj^ is the initial quasi-momentum considering the fi- 
nite size of the BEG, u = f cos^ + 1 j , and 
K {u) is the complete elliptic integral of the first kind as 


K(«) = I 1+ E 




2’^n! 


,2n 


It is obvious that 


t n=l -■ J 

the periods for the two kinds of oscillations are inversely 
proportional to uJx. 

This inverse relation does not change even when we 
use the realistic Bloch bands instead of the idealized 
cosine form. This is confirmed by our numerical results 
with the semi-classical equations (1,2) as shown in 
Fig. 4. Since = 27r x 28Hz in our experiment, we 
have Tp = 17.1ms for Vb = ISE’r and Tpo = 23.9ms for 
Vb = bEr^ which agree very well with the experimental 
results in Fig. 2. In Fig. 3(a), the semi-classical oscil¬ 
lations are plotted as white circles, matching both the 
experimental result and the numerical GPE result. 


V. HIGH ORDER BRAGG REFLECTION 

Bragg reflection (or scattering) at the Bloch band edge 
or center is a fundamental quantum process in periodic 
systems. As we have already seen, it plays a crucial role 
in the oscillations observed in our experiments. We now 
take a closer look at it by recording the absorption im¬ 
ages every 0.1ms. Two series of images are shown in 
Fig. 5, where we see clearly the reflection process be¬ 
tween ±3fi/cL in (a) under IbEr and between ± 2 hkL 
in (b) under bEr. Bragg reflections have been demon¬ 
strated in lower bands with ultracold [1, 33, 34] atomic 
gases, and in high bands with ultracold atoms by using 
a time-dependent optical lattice in Ref. [35] , where only 
reflection results were given, but the reflection processes 
weren’t shown. To the best of our knowledge, we are 
the first to directly observe high order Bragg reflections 
process in high excited bands of optical lattices. 

We use Fig. 5(b) to show in detail what is observed 
in our experiment. The Bragg reflection occurs roughly 
between 5.8ms and 6.1ms. During this period, atoms 
around 2 hkL get scattered to — 2 hkL while atoms around 
—2hkL get scattered to 2hkL. As a result, there are two 
fractions of atoms at 2 hkL' one fraction waits to be re¬ 
flected to — 2 hkL and the other consists of atoms just 
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FIG. 4: (Color online) (a) The periods of the two kinds of 
oscillations for: Tf for Vo = 15F^r; (b) Tfd for 5Er. The 
red lines are fitting functions: (a) Tf = S.Ol/cJcc; (b) Tfd = 
1.38/ujx- 
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scattered from —2hkL. As atoms at ±.2hkL are located 
in different places, these two fractions at 2 hkL are sep¬ 
arated in coordinate space. With a 28ms TOF, these 
two fractions appear as two different peaks. There are 
also two fractions ioi —2hkL- As a result, we observe 
four peaks during the period between 5.8ms and 6.1ms, 
instead of two peaks for other times. 

As shown in Fig. 5(a), the reflection between ±.8hkL 
occurs roughly from 3.4ms to 4.0ms, which is about twice 
longer than the reflection between E2hkL- This shows 
that higher-order Bragg reflection is more difficult to oc¬ 
cur. Fig. 5(c) is the simulation result with the GPE for 
bEr which is consistent with our experimental result. 

When the atoms are not Bragg reflected, they tunnel 
from one band to a neighboring band. This is the well 
known LZ tunneling. In other words, the Bragg reflection 
can be described as a complementary process to the LZ 
tunneling. Between the F and D bands, the tunneling is 
determined by the PpD = [36] with 


7fd = 


^FD 


4h 


6hkLf{x) 


-1 


(5) 


The calculated probability Pfd versus Vb as blue solid 
line is shown in Fig. 6. It is clear that the tunneling prob¬ 
ability decreases with the increasing lattice depth. It is 
close to 1 at Vb = bE^ while close to zero at 15P^. Black 
circles are experimental results for the ratio of atoms tun¬ 


FIG. 5: (Golor online) High order Bragg reflections process 
observed (a) at the edge of the F-band for Vb = 15Pr; (b) at 
the center of the D-band for Vb = 5Pr in the FBZ. (c) is the 
simulation results from GPE for Vb = 5Pr- 


neled to the D-band. There is a very good agreement 
between the experiment and the theory. The band gaps 
Afd between the E and D bands at quasi-momenta be¬ 
ing EhkL versus Vb are shown as red solid line in Eig. 6 
to assist the understanding of this quantum transition. 


VI. CONCLUSION 

In summary, we have loaded a BEC in a one¬ 
dimensional optical lattice non-adiabatically into the G- 
band. The BEC is then observed to tunnel to the E- 
band and undergoes quantum oscillations within the E 
and D bands. The variable force exerted on atoms en¬ 
ables a direct visualization of oscillations in high bands 
in the momentum space. By controlling the Bragg re¬ 
flection at the edges of the E and D bands with optical 
lattice, we have observed three different types of quan¬ 
tum oscillations. At weak lattice strength, oscillations 
between both the E and D bands are observed; at strong 
lattice strength, oscillations only within the E-band are 
observed. At intermediate strength, a superposition of 
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FIG. 6: (Color online) The calculated transition probability 
from the F-band to D-band versus Vo is shown as blue solid 
line. The black circles are the experimental results. The band 
gaps between the F and D bands at the FEZ edge versus the 
lattice depth are also shown as red solid line. 


the above two dynamical oscillations is observed. Fur¬ 
thermore, we has directly demonstrated the high order 
Bragg reflections process in high excited bands of optical 
lattices. 
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